B-Cell and Monocyte Contribution to Systemic Lupus Erythematosus Identified by Cell-Type-Specific Differential Expression Analysis in RNA-Seq Data

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by complex interplay among immune cell types. SLE activity is experimentally assessed by several blood tests, including gene expression profiling of heterogeneous populations of cells in peripheral blood. To better understand the contribution of different cell types in SLE pathogenesis, we applied the two methods in cell-type-specific differential expression analysis, csSAM and DSection, to identify cell-type-specific gene expression differences in heterogeneous gene expression measures obtained using RNA-seq technology. We identified B-cell-, monocyte-, and neutrophil-specific gene expression differences. Immunoglobulin-coding gene expression was altered in B-cells, while a ribosomal signature was prominent in monocytes. On the contrary, genes differentially expressed in the heterogeneous mixture of cells did not show any functional enrichment. Our results identify antigen binding and structural constituents of ribosomes as functions altered by B-cell- and monocyte-specific gene expression differences, respectively. Finally, these results position both csSAM and DSection methods as viable techniques for celltype-specific differential expression analysis, which may help uncover pathogenic, cell-type-specific processes in SLE

Epstein Barr Virus Interleukin 10 Suppresses Anti-inflammatory Phenotype in Human Monocytes

Epstein Barr virus (EBV) is a gamma herpes virus associated with certain malignancies and autoimmune diseases. EBV maintains latency in B cells with occasional reactivation, in part by overcoming the host immune response with viral homologs of several human proteins. EBV interleukin 10 (vIL-10), a lytic phase protein, is a homolog of human IL-10 (hIL-10). The effect of vIL-10 on human monocytes, which are one of the first immune cells to respond to infection, is not known. To understand the role of vIL-10, monocytes from peripheral blood mononuclear cells were stimulated with hIL-10 or vIL-10. Human IL-10 stimulated STAT3 phosphorylation, which is required for suppression of inflammatory responses. However, vIL-10 induced significantly lower phosphorylation of STAT3 compared to hIL-10, and was less efficient in downregulating inflammatory genes. vIL-10 significantly reduced the expression of scavenger receptor CD163 on monocytes, suggesting inhibition of M2 polarization. Furthermore, uptake of apoptotic cells was reduced in vIL-10-stimulated monocytes compared to hIL-10-stimulated monocytes. A neutralizing antibody to IL-10R1 inhibited STAT3 phosphorylation induced by either hIL-10 or vIL-10, suggesting that vIL-10 signals through IL-10R1. Interestingly, vIL-10 suppressed hIL-10-induced STAT3 phosphorylation and inhibited upregulation of suppressors of inflammatory response by hIL-10. We further show that vIL-10 levels were significantly higher in plasma samples from systemic lupus erythematosus (SLE) patients compared to matched unaffected controls. vIL-10 levels did not correlate with hIL-10 levels, but were associated with levels of IgA antibodies to EBV viral capsid antigen, which is an indirect measure of viral reactivation. We propose that the suppression of hIL-10- induced anti-inflammatory genes by vIL-10, together with an increase in inflammatory gene expression, may overcome the anti-inflammatory effects of hIL-10 and exacerbate autoimmune responses in systemic autoimmune diseases

Identification of Unique MicroRNA Signature Associated with Lupus Nephritis

MicroRNAs (miRNA) have emerged as an important new class of modulators of gene expression. In this study we investigated miRNA that are differentially expressed in lupus nephritis. Microarray technology was used to investigate differentially expressed miRNA in peripheral blood mononuclear cells (PBMCs) and Epstein-Barr Virus (EBV)-transformed cell lines obtained from lupus nephritis affected patients and unaffected controls. TaqMan-based stem-loop real-time polymerase chain reaction was used for validation. Microarray analysis of miRNA expressed in both African American (AA) and European American (EA) derived lupus nephritis samples revealed 29 and 50 differentially expressed miRNA, respectively, of 850 tested. There were 18 miRNA that were differentially expressed in both racial groups. When samples from both racial groups and different specimen types were considered, there were 5 primary miRNA that were differentially expressed. We have identified 5 miRNA; hsa-miR-371-5P, hsa-miR-423-5P, hsa-miR-638, hsa-miR-1224-3P and hsa-miR-663 that were differentially expressed in lupus nephritis across different racial groups and all specimen types tested. Hsa-miR-371-5P, hsa-miR-1224-3P and hsa-miR-423-5P, are reported here for the first time to be associated with lupus nephritis. Our work establishes EBV-transformed B cell lines as a useful model for the discovery of miRNA as biomarkers for SLE. Based on these findings, we postulate that these differentially expressed miRNA may be potential novel biomarkers for SLE as well as help elucidate pathogenic mechanisms of lupus nephritis. The investigation of miRNA profiles in SLE may lead to the discovery and development of novel methods to diagnosis, treat and prevent SLE

Internal standard-based analysis of microarray data2—Analysis of functional associations between HVE-genes

In this work we apply the Internal Standard-based analytical approach that we described in an earlier communication and here we demonstrate experimental results on functional associations among the hypervariably-expressed genes (HVE-genes). Our working assumption was that those genetic components, which initiate the disease, involve HVE-genes for which the level of expression is undistinguishable among healthy individuals and individuals with pathology. We show that analysis of the functional associations of the HVE-genes is indeed suitable to revealing disease-specific differences. We show also that another possible exploit of HVE-genes for characterization of pathological alterations is by using multivariate classification methods. This in turn offers important clues on naturally occurring dynamic processes in the organism and is further used for dynamic discrimination of groups of compared samples. We conclude that our approach can uncover principally new collective differences that cannot be discerned by individual gene analysi

Effects of IRF5 Lupus Risk Haplotype on Pathways Predicted to Influence B Cell Functions

Both genetic and environmental interactions affect systemic lupus erythematosus (SLE) development and pathogenesis. One known genetic factor associated with lupus is a haplotype of the interferon regulatory factor 5 (IRF5) gene. Analysis of global gene expression microarray data using gene set enrichment analysis identified multiple interferon- and inflammation-related gene sets significantly overrepresented in cells with the risk haplotype. Pathway analysis using expressed genes from the significant gene sets impacted by the IRF5 risk haplotype confirmed significant correlation with the interferon pathway, Toll-like receptor pathway, and the B-cell receptor pathway. SLE patients with the IRF5 risk haplotype have a heightened interferon signature, even in an unstimulated state (P = 0.011), while patients with the IRF5 protective haplotype have a B cell interferon signature similar to that of controls. These results identify multiple genes in functionally significant pathways which are affected by IRF5 genotype. They also establish the IRF5 risk haplotype as a key determinant of not only the interferon response, but also other B-cell pathways involved in SLE

Tubular cell and keratinocyte single-cell transcriptomics applied to lupus nephritis reveal type I IFN and fibrosis relevant pathways.

The molecular and cellular processes that lead to renal damage and to the heterogeneity of lupus nephritis (LN) are not well understood. We applied single-cell RNA sequencing (scRNA-seq) to renal biopsies from patients with LN and evaluated skin biopsies as a potential source of diagnostic and prognostic markers of renal disease. Type I interferon (IFN)-response signatures in tubular cells and keratinocytes distinguished patients with LN from healthy control subjects. Moreover, a high IFN-response signature and fibrotic signature in tubular cells were each associated with failure to respond to treatment. Analysis of tubular cells from patients with proliferative, membranous and mixed LN indicated pathways relevant to inflammation and fibrosis, which offer insight into their histologic differences. In summary, we applied scRNA-seq to LN to deconstruct its heterogeneity and identify novel targets for personalized approaches to therapy

Neutrophils isolated from systemic lupus erythematosus patients exhibit a distinct functional phenotype

Neutrophil dysregulation, particularly of a low-density subset, is associated with systemic lupus erythematosus (SLE); however, the exact role of normal-density neutrophils in SLE remains unknown. This study compares activation and functional phenotypes of neutrophils from SLE patients and healthy controls to determine potential contributions to SLE pathogenesis. Surface activation markers and release of neutrophil extracellular traps (NETs), granule proteins, and cytokines/chemokines were measured in resting and stimulated neutrophils from SLE patients (n=19) and healthy controls (n=10). Select miRNA and mRNA involved in neutrophil development and function were also measured. Resting SLE neutrophils exhibited fewer activation markers compared to control neutrophils, and activation markers were associated with different plasma cytokines/chemokines in SLE patients compared to healthy controls. However, activation markers increased similarly in SLE and control neutrophils following stimulation with a TLR7/8 agonist, neutrophil growth factors, and bacterial mimic. At the resting state, SLE neutrophils produced significantly more CXCL10 (IP-10), with trends toward other increased cytokines/chemokines. Following stimulation, SLE neutrophils produced fewer NETs and proinflammatory cytokines compared to control neutrophils but more MMP-8. In addition, SLE neutrophils expressed less miR130a, miR132, miR27a, and miR223. In conclusion, SLE neutrophils exhibit distinct functional responses compared to control neutrophils. These functional differences may result from differential gene expression via miRNAs. Furthermore, the differences in functional phenotype of SLE neutrophils suggest that they may contribute to SLE differently dependent on the inflammatory milieu

Immune profiling in Puerto Rican injection drug users with and without HIV-1 infection

Antiretroviral therapy has been effective in suppressing HIV viral load and enabling people living with HIV to experience longer, more conventional lives. However, as people living with HIV are living longer, they are developing aging-related diseases prematurely and are more susceptible to comorbidities that have been linked to chronic inflammation. Coincident with HIV infection and aging, drug abuse has also been independently associated with gut dysbiosis, microbial translocation, and inflammation. Here, we hypothesized that injection drug use would exacerbate HIV-induced immune activation and inflammation, thereby intensifying immune dysfunction. We recruited 50 individuals not using injection drugs (36/50 HIV+) and 47 people who inject drugs (PWID, 12/47 HIV+). All but 3 of the HIV+ subjects were on antiretroviral therapy. Plasma immune profiles were characterized by immunoproteomics, and cellular immunophenotypes were assessed using mass cytometry. The immune profiles of HIV+/PWID−, HIV−/PWID+, and HIV+/PWID+ were each significantly different from controls; however, few differences between these groups were detected, and only 3 inflammatory mediators and 2 immune cell populations demonstrated a combinatorial effect of injection drug use and HIV infection. In conclusion, a comprehensive analysis of inflammatory mediators and cell immunophenotypes revealed remarkably similar patterns of immune dysfunction in HIV-infected individuals and in people who inject drugs with and without HIV-1 infection

Longitudinal patterns and predictors of response to standard-of-care therapy in lupus nephritis: data from the Accelerating Medicines Partnership Lupus Network.

BACKGROUND: Leveraging the Accelerating Medicines Partnership (AMP) Lupus Nephritis (LN) dataset, we evaluated longitudinal patterns, rates, and predictors of response to standard-of-care therapy in patients with lupus nephritis. METHODS: Patients from US academic medical centers with class III, IV, and/or V LN and a baseline urine protein/creatinine (UPCR) ratio ≥ 1.0 (n = 180) were eligible for this analysis. Complete response (CR) required the following: (1) UPCR < 0.5; (2) normal serum creatinine (≤ 1.3 mg/dL) or, if abnormal, ≤ 125% of baseline; and (3) prednisone ≤ 10 mg/day. Partial response (PR) required the following: (1) > 50% reduction in UPCR; (2) normal serum creatinine or, if abnormal, ≤ 125% of baseline; and (3) prednisone dose ≤ 15 mg/day. RESULTS: Response rates to the standard of care at week 52 were CR = 22.2%; PR = 21.7%; non-responder (NR) = 41.7%, and not determined (ND) = 14.4%. Only 8/180 (4.4%) patients had a week 12 CR sustained through week 52. Eighteen (10%) patients attained a week 12 PR or CR and sustained their responses through week 52 and 47 (26.1%) patients achieved sustained PR or CR at weeks 26 and 52. Week 52 CR or PR attainment was associated with baseline UPCR > 3 (ORadj = 3.71 [95%CI = 1.34-10.24]; p = 0.012), > 25% decrease in UPCR from baseline to week 12 (ORadj = 2.61 [95%CI = 1.07-6.41]; p = 0.036), lower chronicity index (ORadj = 1.33 per unit decrease [95%CI = 1.10-1.62]; p = 0.003), and positive anti-dsDNA antibody (ORadj = 2.61 [95%CI = 0.93-7.33]; p = 0.069). CONCLUSIONS: CR and PR rates at week 52 were consistent with the standard-of-care response rates observed in prospective registrational LN trials. Low sustained response rates underscore the need for more efficacious therapies and highlight how critically important it is to understand the molecular pathways associated with response and non-response

A Comprehensive and Universal Method for Assessing the Performance of Differential Gene Expression Analyses

The number of methods for pre-processing and analysis of gene expression data continues to increase, often making it difficult to select the most appropriate approach. We present a simple procedure for comparative estimation of a variety of methods for microarray data pre-processing and analysis. Our approach is based on the use of real microarray data in which controlled fold changes are introduced into 20% of the data to provide a metric for comparison with the unmodified data. The data modifications can be easily applied to raw data measured with any technological platform and retains all the complex structures and statistical characteristics of the real-world data. The power of the method is illustrated by its application to the quantitative comparison of different methods of normalization and analysis of microarray data. Our results demonstrate that the method of controlled modifications of real experimental data provides a simple tool for assessing the performance of data preprocessing and analysis methods